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ABSTRACT:  We  report  a  density  functional  theory  (DFT)  and 
time-dependent  DFT  (TDDFT)  investigation  of  the  thiolated  silver 
nanoclusters  [Ag^SR)^]4-,  Ag14(SR)12(PR'3)8,  Ag31(SG)19/  Ag32- 
(SG)19,  and  Ag15(SG)u,  which  were  synthesized  and  for  which 
one-photon  absorption  (OPA)  characterization  is  available.  Our 
computational  investigation  based  on  careful  examination  of  the 
exchange-correlation  functional  used  in  DFT  geometry  optimization 
and  for  the  linear  optical  properties  predictions  by  TDDFT,  demon¬ 
strated  good  agreement  with  the  measured  linear  absorption 
spectra,  however  dependent  on  the  applied  functional.  Following 
the  benchmarking,  we  evaluated  the  two-photon  absorption  (TPA) 
response  using  TDDFT,  noting  that  accurate  prediction  of  OPA  is 
important  for  suppositions  on  the  spectral  range  for  TPA 
enhancement  because  of  the  sensitivity  to  the  excitation  energies. 

Although  the  TPA  cross-section  results  are  complicated  by  resonance  effects  and  quantifying  TPA  cross  sections  for  these 
systems  is  difficult,  our  results  indicate  that  the  nanoclusters  Ag15  and  Ag31/32  are  likely  to  have  large  TPA  cross  sections.  The 
spherical  symmetry  of  the  Ag44  and  Ag14  nanoclusters  leads  to  applicability  of  superatom  theory,  while  it  is  not  as  useful  for  the 
more  oblate  geometries  of  the  Agls  and  Ag31/32  systems. 


■  INTRODUCTION 

Properties  of  solution-phase  nanoclusters  in  the  size  regime  of 
1—3  nm  have  drawn  significant  interest  because  of  increased 
surface  area  and  changes  in  the  electronic  structure  as  dependent 
on  the  bonding  structure,  among  other  attributes.1  Applications 
ranging  from  enhanced  catalysis,  magnetic  resonant  imaging 
contrast  agents,  energy  conversion,  reduced  toxicity,  and  bio¬ 
logical  sensing  and  imaging  were  investigated.  Furthermore, 
nonlinear  optical  properties  of  thiolated  gold  nanoclusters 
demonstrated  remarkable  nonlinear  optical  enhancements, 
for  example,  the  two-photon  absorption  (TPA)  cross-section 
(<7^)  (imaginary  part  of  second-order  hyperpolarizability)  for  a 
thiolated  Au25(SR)18_1  nanocluster  was  reported  as  large  as 
42  700  GM  at  800  nm,2'3  surpassing  TPA  cross  sections  of 
organic  compounds,  such  as  the  DANS  chromophore,  for  which 
in  dimethyl  sulfoxide  is  191  GM.  Very  large  TPA 
cross  sections  were  more  recently  measured  for  thiolated  gold 
nano  clusters4-0  and  also  predicted  by  calculation.7  Large  static 
first-order  hyperpolarizabilites  for  thiolate-protected  gold  nano¬ 
clusters  were  also  calculated.8  At  the  same  time,  although  optical 
properties  of  thiolated  silver  nanoclusters  were  found  useful  in 
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biological  applications,9,10  a  theoretical  investigation  of  the 
nonlinear  optical  properties  of  these  materials  is  lacking. 

Stable  thiolated  silver  nanoclusters  with  spherical  symme¬ 
try  such  as  [Ag44(SR)30]4-  and  Ag14(SR)12(PR'3)8  have  been 
synthesized  and  characterized  in  the  linear  optical  regime.11-16 
In  the  experimental  work  on  [Ag44(SR)30]4-,  the  organic  group  R 
was  either  benzoic  acid  or  a  fluorinated  phenyl  group,  while  for 
Ag14(SR)12(PR'3)8,  R  =  C6H3F2  and  R'  =  C6H5.  These  materials 
tend  to  follow  the  superatom  number  rule,  often  useful  for  the 
prediction  and  interpretation  of  the  stability  and  properties  of 
metal  nanoclusters.  In  superatom  theory,  the  valence  electrons 
are  assumed  to  behave  as  free  electrons  and,  particularly  if 
the  system  is  approximately  spherical,  will  fill  electron  shells 
analogous  to  those  in  an  atom.  In  a  ligated  system  such  as  the 
gold  or  silver  thiolates,  the  number  of  valence  electrons  ft*  is 
given  by  ft*  =  NvA  —  M  —  z,  where  N  is  the  number  of  metal 
atoms,  vA  is  the  valence  of  the  metal  atom  (assumed  to  be  1  for 
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gold  or  silver)^  M  is  the  number  of  electron-withdrawing  ligands, 
and  z  is  the  net  charge  of  the  molecule.  For  nanoclusters  the 
shell  filling  order  is  1S21 1P61 1D10I2S21F14I2P61G18,  yielding  magic 
numbers  for  stability  of  2,  8,  18,  34,  58,  etc.  (for  some  systems 
20  and  40  may  also  be  part  of  this  sequence). 

Recent  experimental  evidence  has  shown  that  stable  thiolated 
gold  nanoclusters  with  less  symmetry  are  not  well-described  by 
superatom  theory,  and  indeed  a  concern  has  been  raised  by  Joshi 
et  al.17  regarding  the  development  of  a  systematic  approach  to 
the  study  of  ligand  protected  metal  nanoclusters.  Interestingly, 
while  systems  with  lower  symmetry  may  have  less  distinct  one- 
photon  and  two-photon  states,  there  could  be  the  potential  for 
even  larger  TPA  cross  sections,  similar  to  the  larger  TPA  cross 
sections  obtained  from  less  symmetrical  non-centrosymmetric 
donor— acceptor  systems.18 

In  this  work,  to  assess  the  nonlinear  optical  response  of  silver 
nanoclusters,  we  calculated  the  OPA  and  TPA  properties  of 
[Ag44(SR)30]4-  and  Ag14(SR)12(PR'3)8,  and  for  comparison  also 
of  the  Ag31/32(SG)19  and  Agls(SG)n  nanoclusters  by  applying 
DFT  and  TDDFT.  Synthesis  and  optical  characterization  of  the 
silver  glutathione  nanoclusters  Ag32(SG)19  and  Agls(SG)n  were 
recently  reported  by  Ashenfelter  et  al.,19  while  Bertorelle  et  al.20 
reported  experimental  results  for  Ag31(SG)19  and  Agls(SG)u.  In 
these  systems  the  thiol  SG  group  denotes  glutathione,  a 
tripeptide  of  glutamate,  cysteine,  and  glycine,  where  the  sulfur 
atom  on  cysteine  bonds  to  the  silver  atom.  To  mitigate  issues  of 
incorrect  prediction,  careful  examination  of  the  exchange- 
correlation  functional  was  undertaken,  as  used  in  geometry 
optimization,  as  well  as  for  linear  and  nonlinear  optical  properties 
calculations,  including  GGA,  meta-GGA,  hybrid,  and  range- 
separated  hybrid  functionals,  which  we  previously  bench- 
marked.  We  found  that  while  prediction  of  TPA  for  these 
materials  is  a  challenge  due  to  the  importance  of  resonance 
enhancement  and  its  sensitivity  to  excitation  energies,  a  larger 
TPA  cross-section  for  the  less  symmetrical  Ag31/Ag32  and  Ag15 
nano  clusters  was  indicated. 


■  COMPUTATIONAL  METHODS 

The  expression  for  calculating  the  TPA  cross-section  has  been 
previously  derived22-26  (see  references  therein),  given  by 

8tt4 
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where  g  is  a  line  width  function,  \Si0(uvu2)\2  is  the  two-photon 
probability  corresponding  to  a  transition  from  the  ground-state 
(0)  to  a  final  state  (f). 
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where  El  and  E2  are  the  energies  of  the  two  photons  with  unit 
polarization  vectors  ux  and  u2,  respectively.  The  transition  dipole 
moments  are  given  by  the  state  energies  by  E0  and  the  state 
decay  constants  by  IT 

Each  term  in  the  sum-over-states  (SOS)  in  Equation  2  involves 
three  states  (ground-state,  final-state,  and  the  intermediate-state 
over  which  the  sum  is  carried),  and  in  the  essential  three-state 
model,  the  assumption  is  that  a  single  term  dominates  in  the 
SOS.  In  some  cases,  the  TPA  cross-section  can  be  estimated 
accurately  by  this  approximation.  This  is  particularly  true  in 
centrosymmetric  systems  because  the  terms  involving  the  dipole 
moments  are  zero,  and  a  more  limited  number  of  state  transitions 
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are  dipole-allowed.  For  example,  in  a  system  with  C2h  symmetry, 
the  ground  state  has  Ag  symmetry,  the  strongly  allowed  OPA 
states  are  of  Bu  symmetry,  and  the  strong  TPA  states  are  Ag.  In  the 
three-state  approximation  for  the  case  where  the  two  photons 
have  the  same  energy  (E^,),  the  cross-section  is  given  by 
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where  gmax  is  the  maximum  of  the  line  width  function,  and  is 
the  angle  between  the  two  transition  dipole  moment  vectors. 
While  the  three-state  approximation  is  often  not  adequate  for 
quantitative  results,  it  can  be  useful  in  analyzing  the  origin  of  the 
TPA  intensity,  as  we  have  previously  shown.  As  we  have  dis¬ 
cussed  previously,  the  TPA  probability  expressed  in  Equation  2 
can  be  obtained  directly  from  the  residue  of  the  quadratic 
response,  providing  a  more  efficient  method  for  calculating  TPA 
cross  sections.  Because  of  the  size  of  the  clusters  studied  in 
the  paper,  this  is  the  method  we  used  to  calculate  the  TPA 
cross  sections.  This  method  has  the  disadvantage  of  not  having 
the  damping  factor  Ti}  which  can  be  useful  for  mitigating 
unphysically  large  values  near  resonances.  This  damping  factor 
represents  the  line  width  of  the  virtual  state  i,  and  when  used  has 
typically  been  set  to  a  small  value  such  as  0.05  eV.  While  this 
factor  generally  has  little  effect  on  the  TPA  cross-section  for 
organic  systems,  the  higher  density  of  states  in  the  metal  cluster 
systems  makes  near-resonance  conditions  more  likely,  and 
the  inclusion  of  an  accurate  damping  factor  may  dampen  the 
enhancement,  although  it  is  still  expected  to  occur. 

Crystal  structures  for  [Ag44(SR)304-] 11,12  and  [Ag14(SR)12- 
(PR's)*]  15  were  used  as  starting  points  for  geometry  optimiza¬ 
tions  using  DFT.  The  Ag44(SH)304-  system  was  optimized  using 
the  BP86  ,28  exchange-correlation  (X-C)  functional  with  the 
zero-order  regular  approximation  (ZORA)  TZP  basis  set29  in 
ADF.  1,2  The  Ag14(SR)12(PR'3)8  system  was  optimized  with  the 
PW9132  GGA  functional  using  ZORA-TZP29  in  ADF,30'31  as  well 
as  with  the  GGA  PBE3  ■  and  hybrid  PBEO34  (X-C)  functionals 
and  the  LANL2DZ  basis  set35  in  Gaussian  09Z6  The  starting 
structures  used  for  the  Ag31/Ag32  and  Ag15  systems  are  discussed 
in  the  Results  and  Discussion  section.  The  Ag32  system  has  an 
odd  number  of  electrons,  and  the  ground  state  is  a  doublet;  spin- 
unrestricted  DFT  was  used  in  this  case.  DFT  optimizations 
using  the  hybrid  B3LYP27, 37-40  and  PBEO34  functionals  were 
performed  with  the  LANL2DZ  basis  set35  in  Gaussian  09Z6  OPA 
spectra  were  calculated  using  linear  response  TDDFT  in  the 
programs  ADF,26,27  Gaussian  09, 36  GAMESS,41  and  Dalton.42 
Dalton  was  used  to  calculate  TPA  spectra  within  quadratic 
response  TDDFT.  The  LB9443  functional  based  on  an  improved 
model  of  the  potential,  presumably  taking  into  account  the 
asymptotic  behavior,  and  an  improved  version  SAOP,  4  as  well 
as  the  PBE,33  PBEO,34  M06,45  B3LYP,  and  the  hybrid  range- 
separated  CAMB3LYP46  and  mCAMB3LYP47  (X-C)  func¬ 
tionals  were  used  in  the  TDDFT  calculations  with  the  ZORA- 
TZP,29  LANL2DZ,35  and  SDD48’49  basis  sets.  The  polarizable 
continuum  model  (PCM)50  was  used  to  include  solvent  effects 
where  indicated. 


■  RESULTS  AND  DISCUSSION 

A944-  The  structure  and  linear  absorption  spectra  in  the 
[Ag  44  (sr)30]4“  nanocluster  were  reported  by  Desireddy 
et  al., 11,12  and  experimental  and  PBE  spectra  were  also  reported 
by  Yang  et  al.lj  The  experimental  linear  absorption  peaks  are 
at  nearly  the  same  wavelength  whether  the  ligand  (R)  is 
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p-mercaptobenzoic  acid  or  any  of  three  fiuorinated  arylthiols 
(SPhF,  SPhF2,  and  SPhCF3),  and  a  simplified  system  was  used  in 
our  calculations,  with  R  =  H  and  D2h  symmetry  enforced,  as 
shown  in  Figure  1.  The  geometry  was  optimized  with  several 


Figure  1.  Optimized  structure  of  Ag^SR^4  (R  =  H;  see  text  for 
details). 


functionals,  using  the  Desireddy  et  al.11  structure  as  a  guide.  The 
BP86  functional  provided  the  best  agreement  with  experiment; 
thus,  this  optimized  structure  was  used  in  the  TDDFT 
calculations. 

The  structure  of  the  nanocluster  consists  of  a  hollow 
icosahedral  core  of  12  silver  atoms  surrounded  by  20  silver 
atoms  in  the  form  of  a  dodecahedron,  with  each  of  these  20 
atoms  centered  above  one  of  the  20  triangular  faces  of  the  core 
icosahedron,  maintaining  the  Ih  symmetry.  This  structure  is 
surrounded  by  24  sulfur  atoms  in  an  approximate  rhombicu- 
boctahedral  (RCO)  arrangement.  As  this  structure  has  Oh 
symmetry,  this  reduces  the  symmetry  of  the  structure  to  Th. 
The  structure  is  completed  with  six  Ag2S  caps  on  six  faces  of  the 
RCO,  which  maintains  the  approximate  Th  structure.  Alternative 
views  of  the  structure  could  have  Ag2S5  capping  structures  or 
Ag4S5  capping  structures,  where  in  the  latter  view,  two  of  the 
silver  atoms  have  been  extracted  from  the  dodecahedron  for  each 
cap,  and  the  remaining  eight  silver  atoms  in  the  dodecahedron 
form  a  cube  that  encases  the  core  icosahedron,  yielding  an 
alternative  Ag202+  core.  1  While  this  view  of  the  structure  involves 
the  smallest  formal  charges,  the  cubic  core  was  found  to  have  no 
highest  occupied  molecular  orbital— lowest  unoccupied  molec¬ 
ular  orbital  (HOMO— LUMO)  gap.  In  the  case  of  the  Ag2S5 
capping  structures,  however,  a  large  HOMO— LUMO  gap  was 
found  for  the  approximately  spherical  Ag3214+  core.  On  the  basis 
of  this  analysis  of  the  electronic  structure  for  the  two  alternative 
views,  it  was  concluded  that  the  correct  view  of  the  capping 
structures  was  Ag2S5  mounts.11 

In  our  calculations  at  the  BP86/ /TZP-ZORA  level  of  theory 
(see  Computational  Methods  section),  the  silver  atoms  that 
would  be  extracted  from  the  dodecahedron  to  form  Ag4S5  caps 
have  strong,  short  bonds  (2.98  A)  to  the  inner  icosahedral  silver 
atoms,  and  weaker  longer  bonds  (3.33  A)  to  the  outer  Ag  atoms, 
indicating  that  the  preferable  view  is  indeed  that  of  Ag2S5  as 
capping  structures  and  the  dodecahedral  silver  atoms  as 
providing  sturdy  mounts.  Also,  each  of  the  RCO  sulfur  atoms 
is  bonded  to  three  silver  atoms,  with  the  strongest  and  shortest 
bond  (2.56  A)  to  the  outer  Ag,  the  intermediate  length  bond 
(2.65  A)  to  the  mounting  Ag,  and  the  weakest  and  longest  bond 
(2.71  A)  to  the  corner  Ag,  further  supporting  the  view  of  Ag2S5 
caps.  The  outermost  sulfur  atoms  each  have  just  two  bonds,  and 
these  are  to  the  outer  Ag  atoms  and  are  similar  in  length  (2.56  A) 


to  the  other  Ag—S  bonds  in  the  capping  structure.  [Ag44(SR30)]4_ 
crystal  structures  reported  by  Yang  et  al. 1 3  have  the  same  quali¬ 
tative  features  that  support  a  model  with  Ag2S5  caps.  When 
R  =  PhF,  the  bond  lengths  between  the  RCO  sulfur  atoms  and  the 
Ag  atoms  have  the  same  ordering:  2.538  A  (to  outer  Ag  atoms), 
less  than  2.559  A  (to  mount  Ag  atoms),  and  less  than  2.638  A  (to 
corner  Ag  atoms),  although  the  differences  are  not  as  large  as  in 
our  calculations. 

The  excitation  energies  of  the  first  five  peaks  in  the  spectrum, 
including  measured  results  and  those  calculated  with  six  different 
functionals,  as  well  as  the  mean  unsigned  errors  (MUEs)  relative 
to  measurement,  are  listed  in  Table  1.  The  superatom  theory's 


Table  1.  Linear  Absorption  Dataa  for  Ag44(SR)304 


peakl 

peak2 

peak3 

peak4 

peak5 

MUE 

PBE/ZORA:TZP 

1.163 

1.858 

2.263 

2.978 

3.308 

0.245 

LB94/ZORA:TZP 

1.311 

1.851 

2.191 

2.431 

3.066 

0.106 

SAOP/ZORA:TZP 

1.301 

1.911 

2.261 

2.451 

3.096 

0.084 

B3LYP/LANL2DZ 

1.414 

1.914 

2.279 

2.619 

2.894 

0.057 

CAMB3LYP/ 

LANL2DZ 

1.386 

2.076 

2.536 

2.921 

3.161 

0.216 

mCAM/LANL2DZ 

1.373 

1.948 

2.358 

2.688 

3.058 

0.090 

PBE13 

1.300 

1.931 

2.210 

2.480 

2.980 

0.066 

measured13 

1.488 

1.934 

2.317 

2.567 

3.017 

0.026 

measured11 

1.483 

1.921 

2.304 

2.513 

2.972 

0.000 

aIn  electronvolts.  Geometry  optimized  using  BP86/ZORATZP  and 
fwhm  =  0.2  eV. 


so-called  magic  number  n*  for  this  system  is  18,  filling  the 
superatomic  orbitals  1S21 1P6I  ID10.  The  lowest-energy  excitations 
are  dominated  by  transitions  of  ID  electrons  to  the  unfilled  2S 
and  IF  superatomic  orbitals.  This  large  number  of  delocalized 
electrons,  combined  with  the  inner  shell  of  silver  atoms,  which  is 
isolated  from  the  effects  of  the  thiolate  ligands,  allows  for  a 
number  of  low  energy  transitions.  In  the  TDDFT  calculations, 
the  first  peak  in  the  spectrum  falls  in  the  range  of  1.30—1.41  eV, 
just  below  the  measured  value  of  1.48  eV.  The  B3LYP  functional 
performed  best  overall,  with  a  MUE  of  0.057  eV.  Note  that  these 
values  are  lower  than  those  previously  calculated  with  PBE. 
The  OPA  spectrum  calculated  with  B3LYP  is  compared  to  the 
measured  spectrum  in  Figure  2,  and  in  Figure  SI  in  Supporting 
Information,  where  spectra  calculated  using  other  functionals  are 
also  included.  The  B3LYP  functional  was  thus  used  to  calculate 
the  TPA  for  the  first  100  TPA-allowed  states,  which  extends  to 
~2.9  eV.  The  calculated  TPA  spectrum  is  shown  in  Figure  3 
along  with  the  corresponding  OPA. 

The  resonance  enhancement  factor  t  for  TPA  is  defined  by 

t  =  (^  -  0.5 )  .  Thus,  the  calculated  TPA  peak  at  2.75  eV 

\Ef  J 

(corresponding  to  two  901  nm  photons)  is  the  result  of 
resonance  enhancement  from  the  states  that  contribute  to  the 
first  OPA  peak  at  1.33  eV  (934  nm;  see  Figure  3).  Three 
degenerate  states  of  Blu,  B2u,  and  B3u  symmetries  at  an  excitation 
energy  of  1.3276  eV,  each  with  an  oscillator  strength  of  0.0179, 
cause  enhancement  of  two-photon-allowed  states  in  the  energy 
range  of  2.5— 2.8  eV,  particularly  those  of  Ag  symmetry,  which 
give  the  largest  contribution  to  the  TPA  cross-section.  The  trios 
of  degenerate  OPA  states  at  1.8377  and  2.1865  eV  are  expected 
to  cause  significant  TPA  enhancement  near  3.68  and  4.37  eV, 
respectively,  but  our  TPA  calculation  does  not  extend  to  that 
energy,  as  this  spectral  region  is  of  less  interest  for  application. 
In  Figure  S2  we  show  a  plot  of  t  versus  the  transition  energy. 
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Figure  2.  Linear  absorption  spectrum  of  Ag^SR^4  .  Calculated  spectra  for  SR  =  SH  and  a  line  width  of  0.2  eV,  while  the  measured  spectrum  of 
Desireddy  et  al.  had  SR  =  p-mercaptobenzoic  acid. 
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Figure  3.  OPA  and  TPA  spectra  of  Ag^SR^o4-  using  the  B3LYP 
functional,  with  the  TPA  transition  energies  plotted  double  scale  of  the 
OPA  transition  energies,  to  show  the  relationship  between  the  states 
causing  the  resonance  enhancement. 


As  expected,  the  TPA  states  near  the  resonance  transition 
energy  of  2.66  eV  have  a  large  enhancement.  Calculated  TPA 
cross  sections  for  each  TPA-allowed  state  versus  the  resonance 
enhancement  factor  t  for  that  state  (see  Figure  4)  demonstrate 
that  the  resonance  enhancement  is  an  important  factor  in  the 
TPA  cross  section  for  the  silver  nanocluster.  Interestingly, 
although  the  five  states  with  the  largest  TPA  cross  section 
also  have  large  resonant  enhancement  factors,  the  states  with 
the  highest  resonant  enhancement  of  13  277,  a  trio  of  states  at 
2.702  eV,  have  relatively  small  TPA  cross  sections.  A  number  of 
other  states  with  large  resonance  enhancement  factors  have 
insignificant  contributions  to  the  TPA.  Thus,  while  the  resonance 
enhancement  factor  is  important,  it  is  just  one  factor  in  deter¬ 
mining  the  TPA  cross-section,  along  with  the  ground-state  and 
excited-state  transition  dipole  moments. 

Agi4-  The  crystal  structure  and  linear  absorption  spectrum 
of  the  thiolated  silver  nanocluster  Ag14(SR)12(PR/3)8  has  been 
reported  by  Yang  et  al.15  (R  =  C6H3F2  and  R'  =  phenyl),  and 


Gell  et  al.16  performed  some  computational  analysis  of  this 
system.  Additional  studies  were  performed  by  Muniz-Miranda 
et  al.51  and  Goh  et  al.,52  which  have  mentioned  that  this  structure 
satisfies  the  superatom  nanocluster  criteria.  The  optimized  struc¬ 
tures  are  shown  in  Figure  5. 

The  core  of  the  Ag14(SR)12(PR/3)8  nanocluster  structure 
is  comprised  of  six  silver  atoms  located  at  the  vertices  of  an 
icosahedron.  The  remaining  eight  silver  atoms  are  located  at  the 
corners  of  a  cube  that  surrounds  the  icosahedron.  The  12  sulfur 
atoms  are  located  near  the  center  of  each  of  the  12  edges  of  this 
cube  but  are  distorted  so  as  to  bond  with  only  one  of  the  core 
icosahedral  silver  atoms,  while  bonding  with  two  of  the  corner 
atoms.  The  sulfur  atom's  bond  to  the  core  silver  atom  is  the 
shorter,  stronger  bond,  with  an  average  bond  length  of  2.54  A  in 
the  crystal  structure,  compared  to  2.65  A  for  the  outer  Ag— S 
bonds.  Table  SI  in  Supporting  Information  lists  the  bond  lengths 
for  this  structure  with  different  ligands,  including  the  reported 
crystal  structure,15  the  PBE  optimized  structure  of  Gell,  at  al.,16 
and  our  optimized  structures  using  the  PBE,  PBE0,  and  PW91 X-C 
functionals.  While  the  lowest  structural  errors  were  obtained 
when  the  ligands  were  replaced  with  methyl  groups  and  either  the 
PBE/VASP  or  the  PW91/IMCP-SR1  level  of  theory  was  used,  the 
structures  with  just  methyl  groups  did  not  yield  accurate  spectra  in 
the  TDD  FT  calculations  (see  Table  2). 

The  measured  linear  absorption  spectrum,15  shown  in  Figure  6 
and  listed  in  Table  2,  has  a  small  peak  at  2.34  eV  (530  nm)  and  a 
sharper,  more  intense  peak  at  3.37  eV  (368  nm).  The  nano¬ 
cluster  s  smaller  size  compared  to  [Ag44(SR30)]4_  causes  more 
quantum  confinement,  leading  to  larger  transition  energies  and  a 
more  blue-shifted  spectrum.  When  methyl  groups  are  used  as 
ligands,  the  calculated  spectrum  using  TDDFT  is  too  blue- 
shifted,  but  when  the  C6H4F  group  is  used  as  the  sulfur  ligand, 
good  agreement  with  experiment  is  obtained  using  the  LB94 
functional,  and  fair  agreement  is  obtained  using  the  SAOP 
functional.  The  results  using  the  B3LYP,  PBE0,  CAMB3LYP, 
and  M06  X-C  functionals  are  still  too  blue-shifted.  Note  that 
the  superatom  magic  number  for  Ag14(SR)12(PR'3)8  is  2, 
corresponding  to  the  filling  of  the  IS  superatom  orbital. 
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Figure  4.  Calculated  TPA  cross-section  of  Ag44(SR)304  for  each  TPA-allowed  state  vs  its  resonance  enhancement  factor. 


Figure  5.  Optimized  structures  for  Ag14(SCH3  )12(P(CH3)3)8  and 
Agi4(SC6H5F)12(P(CH3)3)8  (see  text  for  details). 


The  TPA  for  this  system  was  calculated  with  methyl  ligands 
using  the  B3LYP  and  CAMB3LYP  functionals.  When  the  B3LYP 
X-C  functional  was  used,  a  peak  TPA  of  ~40  GM  was  obtained  at 


a  transition  energy  of  ~4.1  eV  (corresponding  to  two  605  nm 
photons).  When  the  CAMB3LYP  functional  is  used,  the  calcu¬ 
lated  peak  intensity  is  only  ~20  GM,  and  it  is  blue-shifted  by 
~0.7  eV  (not  shown).  A  blue  shift  in  the  excitation  energy  when 
using  the  range-separated  hybrid  functional  CAMB3LYP  versus 
hybrid  B3LYP  has  been  previously  explained.  The  TPA  and 
OPA  spectra  calculated  with  the  CAMB3LYP  functional  are 
plotted  in  Figure  7.  The  calculated  TPA  for  this  system  is  less  (by 
a  factor  of  ~400)  than  was  calculated  for  the  nanocluster  with  44 
silver  atoms.  The  strong  TPA  is  in  the  range  of  4—5  eV,  but  since 
there  are  no  OPA  states  in  the  2.0— 2.5  eV  range  (the  OPA  states 
are  above  3.0  eV),  the  resonance  enhancement  factor  is  small, 
and  the  TPA  cross-section  is  small.  Note  that  a  large  resonance 
enhancement  occurs  when  the  excitation  energy  of  an  OPA  state 
is  close  to  half  that  of  a  TPA  state.  Because  all  the  OPA  states 
for  the  Ag14  system  are  significantly  higher  in  energy  than 
half  the  TPA  energy  over  the  range  we  are  investigating,  all  the 
enhancement  factors  in  the  TPA  range  of  interest  are  less 
than  20. 


Table  2.  Linear  Absorption  Spectra  for  Ag14(SR)12(PR,3)8  (eV) 


R 

R' 

geometry 

TDDFT 

MUE 

errl 

err2 

pointl 

point2 

ch3 

ch3 

PW91 

LB94/ZORATZP 

0.26 

0.41 

-0.11 

2.75 

3.26 

ch3 

ch3 

PW91 

SAOP/ZORATZP 

0.69 

0.62 

0.76 

2.96 

4.13 

ch3 

ch3 

PBE 

PBE/ZORATZP 

0.72 

0.87 

0.57 

3.20 

3.94 

ch3 

ch3 

PBE 

LB94/ZORATZP 

0.62 

0.74 

0.50 

3.07 

3.86 

ch3 

ch3 

PBE 

SAOP/ZORATZP 

0.21 

0.41 

-0.01 

2.75 

3.36 

ch3 

ch3 

PBE 

CAMB3LYP/SDD 

1.48 

1.26 

1.70 

3.60 

5.07 

c6h4f 

ch3 

PBE 

LB94/ZORATZP 

0.03 

0.00 

-0.07 

2.34 

3.30 

C6H4F 

ch3 

PBE 

SAOP/ZORATZP 

0.12 

0.16 

0.08 

2.49 

3.44 

c6h4f 

ch3 

PBE 

B3LYP/SDD 

0.64 

0.73 

0.54 

3.07 

3.91 

c6h4f 

ch3 

PBE 

CAMB3LYP/SDD 

1.48 

1.26 

1.70 

3.60 

5.07 

c6h4f 

ch3 

PBE 

M06/LANL2DZ 

0.72 

0.74 

0.71 

3.08 

4.08 

c6h4f 

ch3 

PBE0 

PBE0/LANL2DZ 

1.04 

0.90 

1.17 

3.24 

4.54 

C6H4F 

ch3 

PBE0 

LB94/ZORATZP 

0.06 

0.02 

-0.10 

2.36 

3.27 

c6h4f 

ch3 

PBE0 

SAOP/ZORATZP 

0.10 

0.17 

-0.04 

2.50 

3.32 

c6h4f 

ch3 

PBE0 

B3LYP/LANL2DZ 

0.83 

0.78 

0.87 

3.12 

4.24 

c6h4f 

ch3 

PBE0 

CAMB3LYP/SDD 

1.58 

1.33 

1.82 

3.67 

5.19 

c6h4f 

ch3 

PBE0 

M06/LANL2DZ 

0.80 

0.74 

0.86 

3.08 

4.23 

C6H3F2 

phenyl 

measured15 

2.34 

3.37 
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Figure  6.  Linear  absorption  of  Ag14(SR)12(PR'3)8.  For  calculated  spectra,  R=  C6H4F,  R'  =  CH3,  the  geometry  was  optimized  with  the  PBE  functional, 
and  the  line  width  is  0.5  eV.  For  the  measured  data  from  Yang,  et  al., 15  R  =  C6H3F2  and  R'  =  phenol. 


CAMB3LYP:  - OPA  — TPA 


OPA  transition  energy  (eV) 

TPA  transition  energy  (eV) 

Figure  7.  OPA  and  TPA  spectra  of  Ag14(SR)12(PR'3)8,  with  the  TPA  transition  energies  plotted  double  scale  of  the  OPA  transition  energies. 
R  =  R'  =  methyl. 


Ag3i/Ag32  and  Ag15.  Synthesis  and  optical  characterization 
of  the  Ag32(SG)19,  Ag31(SG)19,  and  Agls(SG)n  silver  glutathione 
nanoclusters  have  been  reported.  '  A  structure  for  Ag15(SG)n 
consisting  of  an  Ag8  core  protected  by  two  LI  ligands,  one  L2 
ligand,  and  one  L3  ligand  (LI,  L2,  and  L3,  have  structures  of 
-S(G)-Ag-S(G)-  — S(G)— Ag— S(G)— Ag— S(G)— ,  and 

— S(G)— Ag— S(G)— Ag— S(G)— Ag— S(G)— ,  respectively)  was 
proposed.20  A  structure  was  also  proposed  for  Ag31(SG)19 
consisting  of  an  Ag21  core  protected  by  eight  LI  ligands  and  one 
L2  ligand.  We  propose  a  similar  structure  for  Ag32(SG)19,  where 
two  of  the  LI  ligands  are  connected  by  the  additional  Ag  atom  to 
form  an  L3  ligand,  resulting  in  a  structure  with  six  LI  ligands,  one 
L2  ligand,  and  one  L3  ligand.  The  geometries  for  these  systems, 
with  SCH3  replacing  SG,  were  optimized  using  the  B3LYP 
functional  and  LANL2DZ  basis  set  (see  Computational  methods 
section).  Both  gas-phase  and  aqueous  (using  PCM  to  include  the 


solvent  effect)  geometries  were  obtained,  and  the  solution-phase 
geometries  are  shown  in  Figure  8.  The  coordinates  for  these 
aqueous  systems  optimized  with  the  PBEO  functional  are  listed  in 
Table  S2. 

Figure  9  and  Tables  3  and  4  show  the  measured  and  calculated 
linear  absorption  spectra  for  the  Ag31/Ag32  systems.  The  mea¬ 
sured  linear  absorption  spectrum  of  the  material  that  Bertorelle 
et  al.21  report  to  be  Ag31(SG)19  has  a  sharp  peak  at  2.52  eV,  while 
a  measured  peak  centered  at  2.59  eV  is  shown  in  spectra  by 
Kumar  et  al.,53  as  well  as  by  Ashenfelter  et  al.,19  labeled  “band  6”, 
and  identified  as  the  nanocluster  Ag32(SG)19.  In  the  TDDFT 
calculations  of  Bertorelle  et  al.,20  where  the  SG  ligands  were 
replaced  by  SCH3,  they  obtained  a  peak  at  2.42  eV,  using 
CAMB3LYP  in  the  gas-phase.  Note  that  neutral  Ag31(SG)19  and 
Ag32(SG)19  have  free  electron  counts  of  12  and  13,  respectively, 
which  are  not  magic  numbers.  If  Ag31(SG)19  had  a  charge  of 
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Figure  8.  Solution-phase  optimized  structures  for  Ag15(SCH3)n, 
Ag3i(SCH3)i9;  and  Ag32(SCH3)19. 


either  +4  or  —6  (+5  or  —5  for  Ag32(SG)19),  it  would  have  a 
“superatom”  free  electron  count  of  8  or  18,  which  are  both  magic 
numbers.  However,  attempts  to  optimize  Ag31(SG)19  with 
charges  +4  and  —6  were  not  successful,  and  a  minimum  energy 
structure  has  been  found  only  for  the  neutral  system.  The 
Ag32(SG)19  structure  was  also  optimized  with  zero  charge, 
and  since  this  system  has  an  odd  number  of  electrons,  spin- 
unrestricted  DFT  was  used.  Because  these  structures  do  not  have 
the  spherical  symmetry  of  the  Ag44  and  Ag14  systems,  the 
superatom  magic  number  theory  may  be  less  applicable. 

As  can  be  seen  in  Figure  9,  when  the  linear  absorption 
spectrum  is  calculated  using  the  CAMB3LYP  functional  and 
PCM  solvation  model,  excellent  agreement  is  obtained  with  this 
measured  peak  for  both  Ag31(SCH3)19  (peak  at  2.52  eV)  and 
Ag32(SCH3)19  (peak  at  2.49  eV),  assuming  a  peak  full-width  at 
half-maximum  (fwhm)  of  0.2  eV  for  all  calculated  states  below 
2.7  eV.  For  higher  energy  states,  a  fwhm  of  0.5  eV  was  used  to 
obtain  better  agreement  with  experiment  in  this  part  of  the 
spectrum.  As  can  be  seen  by  the  results  listed  in  Tables  3  and  4, 
good  agreement  is  obtained  whenever  the  CAMB3LYP  func¬ 
tional  is  used,  whether  or  not  the  PCM  solvation  model  is 
included,  while  when  the  B3LYP,  PBEO,  or  mCAMB3LYP 
functionals  were  used,  the  peak  is  significantly  red-shifted. 
While  the  good  agreement  of  the  calculated  results  with 
measurements  gives  credence  to  the  structures  and  methods 
used,  the  similarity  of  the  absorption  spectra  for  Ag31(SG)19 
and  Ag32(SG)19  in  both  the  measured  and  calculated  results 
means  they  cannot  be  used  to  distinguish  between  the  two 
structures. 


Table  3.  Linear  Absorption  Spectra  for  Ag31(SR)19  (eV) 


geometry 

TDDFT 

error 

peak 

PBE0/ZORA:TZP 

CAMB3LYP/LANL2DZ 

-0.115 

2.400 

B3LYP/LANL2DZ 

PBE0/LANL2DZ 

-0.259 

2.256 

B3LYP/LANL2DZ 

B3LYP/LANL2DZ 

-0.345 

2.170 

B3LYP/LANL2DZ 

CAMB3LYP/LANL2DZ 

0.072 

2.587 

B3LYP/LANL2DZ 

CAMB3LYP/ 

LANL2DZ(PCM) 

0.026 

2.541 

B3LYP/LANL2DZ(PCM) 

PBE0/LANL2DZ(PCM) 

-0.250 

2.265 

B3LYP/LANL2DZ(PCM) 

B3LYP/LANL2DZ(PCM) 

-0.261 

2.255 

B3LYP/LANL2DZ(PCM) 

CAMB3LYP/ 

LANL2DZ(PCM) 

0.008 

2.523 

PBEO /LANL2DZ 

PBE0/LANL2DZ 

-0.274 

2.241 

PBE0/LANL2DZ 

B3LYP/LANL2DZ 

-0.345 

2.170 

PBE0/LANL2DZ 

CAMB3LYP/LANL2DZ 

0.056 

2.571 

PBE0/LANL2DZ(PCM) 

PBE0/LANL2DZ  (  P  CM) 

-0.244 

2.271 

PBE0/LANL2DZ(PCM) 

B3LYP/LANL2DZ(PCM) 

-0.291 

2.224 

PBE0/LANL2DZ(PCM) 

CAMB3LYP/ 

LANL2DZ(PCM) 

-0.028 

2.487 

PBE0/LANL2DZ(PCM) 

mCAMB3LYP/ 

LANL2DZ(PCM) 

-0.199 

2.316 

measured20 

2.515 

The  TPA  for  Ag31(SCH3)19  was  calculated  at  the  CAM- 
B3LYP//PBE0  level  of  theory  and  plotted  in  Figure  10,  along 
with  the  OPA  at  the  same  level  of  theory.  The  TPA  transition 
energies  are  plotted  double  scale  to  aid  in  noticing  resonance 
effects.  As  the  oscillator  strengths  for  the  first  three  excited  states 
are  small  (S1  =  0.89  eV,f=  0.001;  S2  =  1.31  eV,f=  0.001;  S3  = 
1.50  eV,/=  0.004),  the  correspondence  is  not  easily  observed  in 
this  plot.  However,  in  Figure  S3,  the  TPA  from  each  final  state  is 
plotted  versus  the  resonance  enhancement  factor  t  using  each  of 
the  first  three  OPA  states  as  the  intermediate  state.  This  plot 
shows  that  the  TPA  peak  near  2.62  eV  of  668  GM  is  the  result  of 
two  final  states  with  contributions  of  458  GM  and  187  GM  that 
benefit  from  large  resonance  enhancement  factors  of  1  000  000 
and  400  000  due  to  S2.  One  can  also  observe  that  the  small 
TPA  peak  of  17  GM  near  1.9  eV  is  the  result  of  resonance 
enhancement  with  Sx  producing  a  factor  of  t  =  31  000.  Currently 
we  were  able  to  perform  TPA  calculations  on  this  system  for  the 
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Figure  9.  OPA  spectra  for  Ag31/Ag32  nanoclusters. 
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Table  4.  Linear  Absorption  Spectra  for  Ag32(SR)19  (eV) 


geometry 

B3LYP/LANL2DZ(PCM) 

PBE0/LANL2DZ 

PBE0/LANL2DZ 

PBE0/LANL2DZ 

PBE0/LANL2DZ 

PBE0/LANL2DZ(PCM) 

PBE0/LANL2DZ(PCM) 

PBE0/LANL2DZ(PCM) 


TDDFT  error  peak 

CAMB3LYP/LANL2DZ(PCM)  -0.095  2.488 

PBE0/LANL2DZ  -0.451  2.132 

B3LYP/LANL2DZ  -0.546  2.037 

CAMB3LYP/LANL2DZ  -0.097  2.486 

mCAMB3LYP/LANL2DZ  -0.338  2.245 

PBE0/LANL2DZ ( P CM)  -0.259  2.324 

B 3LYP / LANL2DZ ( P CM)  -0.329  2.254 

mCAMB3LYP/LANL2DZ(PCM)  -0.182  2.401 

measured19  2.583 


first  30  states,  extending  to  3.1  eV,  but  we  expect  that  much  larger 
TPA  cross  sections  will  be  obtained  at  transition  energies  greater 
than  4.0  eV  due  to  resonance  enhancement  with  the  stronger 
OPA  states  in  the  range  of  2.3— 2.6  eV. 

The  measured  and  calculated  linear  absorption  spectrum 
for  Agls(SR)u  is  given  in  Table  5.  Figure  11  compares  spectra 
computed  with  the  PBE0  functional  to  the  measured  data,  while 
Figure  S4  includes  calculated  data  using  additional  function¬ 
als.  The  experimental  linear  absorption  spectra  reported  for 
Ag15(SG)u^20  are  distinguished  by  a  steep  rise  between  2.2  and 
2.5  eV  followed  by  a  flat  plateau  from  2.6  to  3.2  eV.  The 
theoretical  spectrum  previously  reported  by  Bertorelle  et  al.20  has 
a  narrow  peak  at  2.17  eV,  another  peak  at  2.85  eV,  and  a  more 
intense  narrow  band  at  3.02  eV,  and  correspondence  with 
measurement  is  not  apparent.  According  to  superatom  theory, 
with  a  neutral  charge  this  system  has  four  free  electrons,  which  is 
not  a  magic  number  for  stability.  In  our  calculations,  both  the 
neutral  system  and  the  cation  with  a  charge  of  +2,  a  system  for 
which  the  number  of  free  electrons  is  the  magic  number  2,  were 
optimized.  However,  spectra  calculated  from  the  cation  are  not  in 
as  good  agreement  with  measurement,  and  we  focus  on  the 
neutral  system  in  the  following.  As  with  the  Ag31  nano  cluster, 
Agls  is  not  spherically  symmetrical,  and  superatom  theory  may 
not  be  applicable.  Our  calculated  spectra  in  Figure  1 1  use  a  line 
width  of  0.5  eV  and  include  both  gas-phase  and  aqueous  results, 
using  the  PCM  solvation  model. 


The  most  distinguishing  feature  in  the  absorption  spectra 
for  Ag15(SG)u19,20  is  the  steep  rise  near  2  eV.  To  quantify  this 
feature,  we  analyzed  the  first  derivative  of  the  absorption 
intensity  with  respect  to  the  excitation  energy,  and  the  maximum 
in  this  derivative  indicates  the  energy  where  the  rise  is  the 
steepest.  This  normalized  derivative  (ND)  is  also  listed  in 
Table  5.  The  results  demonstrate  that,  while  several  calculations 
both  with  and  without  PCM  show  good  agreement  with  experi¬ 
ment  (i.e.,  absolute  error  less  than  0.100  eV),  when  PCM  is  not 
used  the  rise  is  not  steep  enough,  so  the  shape  of  the  plot  is  not  in 
such  good  agreement  with  experiment  (see  Figure  ll).  Overall, 
the  solvated  results  are  in  better  agreement  with  experiment, 
particularly  when  the  PBE0  functional  is  used  for  TDDFT 
calculations.  Fair  results  were  obtained  with  the  B3LYP  and 
mCAMB3LYP  functionals,  while  the  solvated  CAMB3LYP 
results  are  not  in  as  good  agreement  with  experiment  due  to  a 
0.2  eV  blue  shift  of  the  strong  excited  state  relative  to  the 
mCAMB3LYP  calculation. 

Figure  12  shows  the  calculated  TPA  spectra  using  the 
CAMB3LYP  functional  and  two  different  geometries  (from 
B3LYP  and  PBE0  optimizations).  Despite  the  similarities  in  the 
two  functionals  as  well  as  in  the  resulting  structures,  the 
sensitivity  of  the  excitation  energies  to  structure  and  of  the  TPA 
cross  sections  to  resonance  effects  results  in  different  TPA 
spectra.  Calculated  OPA  spectra  are  also  shown  for  these  two 
levels  of  theory,  plotted  at  half-scale  to  illustrate  how  the  single 
photon  resonances  contribute  to  the  calculated  TPA.  When 
the  B3LYP  geometry  is  used,  the  first  excited  state  at  1.85  eV  is  in 
resonance  with  several  states  near  3.7  eV,  resulting  in  an 
extremely  large  calculated  TPA  of  over  60  000  GM.  When  the 
PBE0  functional  is  used,  the  first  excited  state  is  at  1.27  eV, 
but  since  no  calculated  states  are  particularly  close  in  energy  to 
2.54  eV,  there  is  not  a  large  TPA  resonance  effect  from  this  state. 
However,  the  second  excited  state  near  2.11  eV  is  in  resonance 
with  several  states  near  4.22  eV,  and  therefore  this  calculation 
produces  a  large  TPA  at  this  energy.  This  is  further  illustrated  in 
Figure  S5,  where  the  TPA  cross-section  for  each  final  state  using 
the  CAMB3LYP//PBE0  level  of  theory  is  plotted  versus  the 
resonance  enhancement  factor  t,  which  has  been  calculated  using 
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Figure  10.  OPA  and  TPA  spectra  of  Ag31(SR)19,  with  the  TPA  transition  energies  plotted  double  scale  of  the  OPA  transition  energies.  The  line  width  is 
0.2  eV  for  OPA  and  0.3  eV  for  TPA 
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Table  5.  Linear  Absorption  Spectraa  for  Agls(SR)n 


geometry 

TDDFT 

TE 

error 

ND 

B3LYP/LANL2DZ 

SAOP/ZORA:TZP 

2.077 

-0.262 

2.56 

B3LYP/LANL2DZ 

LB94/ZORA:TZP 

1.973 

-0.367 

2.67 

B3LYP/LANL2DZ 

PBE0/LANL2DZ 

2.298 

-0.042 

3.62 

B3LYP/LANL2DZ 

B3LYP/LANL2DZ 

2.207 

-0.132 

3.41 

B3LYP/LANL2DZ 

CAMB3LYP /LANL2DZ 

2.657 

0.318 

3.63 

B3LYP/LANL2DZ 

CAMB3LYP/SDD 

2.606 

0.266 

3.10 

B3LYP/LANL2DZ 

CAMB3LYP  /  SDD 

2.576 

0.236 

3.29 

B3LYP/LANL2DZ 

mCAMB3LYP/LANL2DZ 

2.382 

0.043 

3.72 

PBE0/LANL2DZ 

SAOP/ZORA:TZP 

2.010 

-0.329 

2.74 

PBEO /LANL2DZ 

LB94/ZORA:TZP 

1.982 

-0.357 

2.64 

PBE0/LANL2DZ 

PBE0/LANL2DZ 

2.391 

0.052 

2.39 

PBEO /LANL2DZ 

B3LYP/LANL2DZ 

2.306 

-0.033 

2.57 

PBE0/LANL2DZ 

CAMB3LYP/LANL2DZ 

2.455 

0.116 

2.67 

PBE0/LANL2DZ 

CAMB3LYP/ SDD 

2.340 

0.000 

2.83 

PBEO /LANL2DZ 

mCAMB3LYP /LANL2DZ 

2.450 

0.111 

2.08 

B 3LYP/LANL2D Z ( P CM)  PBE0/LANL2DZ(PCM) 

2.317 

-0.022 

3.87 

B3LYP/LANL2DZ(PCM)  B3LYP/LANL2DZ(PCM) 

2.212 

-0.127 

4.17 

B  3LYP/LANL2D  Z  (  P  CM)  CAMB3LYP/ 

LANL2DZ(PCM) 

2.653 

0.313 

3.90 

B 3LYP/LANL2D Z ( P CM)  mCAMB3LYP/ 

LANL2DZ(PCM) 

2.404 

0.065 

4.17 

PBE0/LANL2DZ(PCM) 

PBE0/LANL2DZ(PCM) 

2.370 

0.030 

3.96 

PBE0/LANL2DZ(PCM) 

B3LYP/LANL2DZ(PCM) 

2.274 

-0.065 

3.87 

PBE0/LANL2DZ(PCM) 

CAMB3LYP/ 

LANL2DZ(PCM) 

2.717 

0.378 

3.95 

PBE0/LANL2DZ(PCM) 

CAMB3LYP/ SDD  (PCM) 

2.440 

0.101 

3.03 

PBE0/LANL2DZ(PCM) 

mCAMB3LYP/ 

LANL2DZ(PCM) 

2.464 

0.125 

3.98 

PBE0/LANL2DZ(PCM) 

mCAMB3LYP/SDD  (PCM) 

2.279 

-0.060 

2.44 

measured(Bertorelle  et  al.)20 

2.458 

0.119 

4.92 

measured(Ashenfelter  et  al.)19 

2.339 

0.000 

5.90 

"Transition  energy  (TE  in  units  of  eV)  where  the  derivative  of  the 
intensity  with  respect  to  the  TE  is  a  maximum,  as  well  as  the  error 
(eV)  compared  to  measurement.  The  normalized  derivative  (ND)  of 
the  intensity  has  units  of  (eV)-1. 


each  of  these  first  two  excited  states  as  the  intermediate  state. 
For  the  final  states  contributing  to  the  very  large  TPA  cross 


section  near  4.2  eV,  the  correlation  is  clearly  with  state  S2,  as  was 
expected  from  Figure  12.  We  emphasize  that  the  level  of  theory 
applied  is  important  for  these  predictions. 

Since  the  calculated  spectra  for  Ag31(SCH3)19  using  the 
CAMB3LYP  functional  are  more  accurate  than  the  results  with 
any  other  functional,  our  discussion  of  the  solvation  effects  on 
this  system  will  be  limited  to  CAMB3LYP.  In  the  calculations  on 
Ag3i(SCFl3)19,  the  effect  of  inclusion  of  the  PCM  solvation 
model  on  the  peak  absorption  energy  is  small,  with  a  red  shift  of 
~0.08  eV  due  to  the  solvent.  The  effect  of  using  the  PCM  for  the 
geometry  optimization  is  small,  with  a  small  average  decrease 
in  the  Ag— Ag  bond  lengths,  and  a  small  average  increase  in  the 
Ag— S  bond  lengths.  From  the  insight  gained  from  the  TDDFT 
calculation  with  PCM  for  the  gas-phase  geometry,  we  can 
attribute  ~30%  of  the  red  shift  to  the  structure  change  and  ~70% 
to  the  use  of  the  PCM  solvation  model  in  the  TDDFT  response 
calculation.  The  effect  of  the  solvation  model  on  the  intensity  is 
more  significant,  with  the  PCM  calculation  yielding  a  peak 
extinction  coefficient  about  twice  that  obtained  in  the  gas-phase 
calculation. 

The  effect  of  the  solvent  model  inclusion  on  the  absorption 
spectra  of  Ag15(SCH3)u  seems  more  significant  since  the  shape 
of  the  measured  absorption  spectra  is  much  better  emulated 
when  the  solvation  model  is  used.  The  best  functional  for  this 
system  is  not  as  obvious,  but  careful  analysis  favors  the  PBEO 
functional.  The  steep  rise  of  the  PCM  spectra  compared  to 
the  gas-phase  data  can  be  mainly  attributed  to  the  much  larger 
oscillator  strength  (J=  0.3 1)  of  the  most  intense  absorption  state 
at  2.55  eV,  compared  to  the  most  intense  state  in  the  gas-phase 
calculation  (2.51  eV,/=  0.09).  Thus,  for  both  systems,  the  biggest 
effect  of  the  solvation  model  is  the  increase  in  intensity  for  the 
most  prominent  absorption  state. 

■  CONCLUSIONS 

Prediction  of  OPA  spectra  in  comparison  with  available 
experimental  data  for  silver  nano  clusters  has  shown  that,  for 
the  OPA  of  [Ag44(SR)30]4_,  the  B3LYP  functional  performed  the 
best,  with  acceptable  performance  from  SAOP  and  mCAM, 
while  for  Ag14(SR)12(PR'3)8,  good  agreement  was  obtained  only 


-  -  PBE0//PBE0  - PBE0(PCM)//PBE0(PCM) 

measured(Ref.  20)  «=measured(Ref.  19) 


Figure  11.  Linear  absorption  of  Ag15(SR)n.  Our  calculated  spectra  used  R  =  CH3  and  a  line  width  of  0.5  eV.  The  measured  results19'20  had 
SR  =  glutathione. 
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■  OPA:CAMB3LYP//PBEO 

■  TPA:CAMB3LYP//PBE0 
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Figure  12.  OPA  and  TPA  of  Ag15(SR)n,  with  the  TPA  transition  energies  plotted  double  scale  of  the  OPA  transition  energies.  Line  width  is  0.2  eV  for 
OPA  and  0.3  eV  for  TPA. 


with  the  LB94  and  SAOP  functionals.  For  the  Ag31/Ag32  systems, 
excellent  results  were  obtained  with  the  CAMB3LYP  functional, 
while  for  the  Agls  system,  good  results  were  obtained  with  the 
PBE0  functional  and  fair  results  were  obtained  with  B3LYP  and 
mCAMB3LYP.  Our  results  indicate  that  using  PBE  for  the 
prediction  of  the  OPA  spectra  of  silver  clusters  may  lead  to  erro¬ 
neous  results.  Note  that  structure  prediction  is  still  a  challenge  for 
silver  nanoclusters,  as  recently  pointed  out  by  Wickramasinghe 
et  al.54 

Indeed,  we  note  that  accurate  prediction  of  OPA  is  important 
for  suppositions  on  the  spectral  range  for  TPA  enhancement 
because  of  the  sensitivity  to  the  excitation  energies.  The  OPA 
spectrum  of  the  Ag44  system  is  largely  independent  of  the  specific 
thiolate  ligand,  likely  due  to  the  isolation  of  the  inner  icosahedral 
shell  of  silver  atoms  from  the  thiolate  groups.  We  found  that  the 
Ag44  and  Ag14  systems  have  significant  spherical  symmetry  and 
follow  the  superatom  rule  for  stability  and  even  for  optical 
properties.  Calculations  on  the  Ag44  and  Agls  systems  predict 
large  TPA,  while  the  Ag14  system  was  found  to  have  a  small  TPA 
over  the  range  calculated  but  may  have  large  TPA  at  higher 
energies  where  resonant  enhancement  may  play  a  bigger  role. 
The  Ag31  system,  which  is  predicted  to  have  a  modest  TPA  cross 
section  of  nearly  700  GM  at  a  transition  energy  near  2.6  eV,  likely 
has  a  much  larger  TPA  cross  section  at  a  transition  energy  closer 
to  5  eV,  due  to  resonance  enhancement.  Finally,  it  is  interesting 
to  note  that,  while  prediction  of  TPA  for  these  materials  is  a 
challenge  due  to  the  importance  of  resonance  enhancement  and 
its  sensitivity  to  excitation  energies,  our  calculations  indicate  large 
TPA  for  the  less  symmetrical  Ag31/Ag32  and  Agls  nanoclusters. 
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